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^ ■ ABSTRACT 
> 

We present sensitive phase-referenced VLBI results on the radio continuum 
emission from the z = 4.4 QSO BRI 1335-0417. The observations were carried 
out at 1.4 GHz using the High Sensitivity Array (HSA). Our sensitive VLBI 
image at 189 x 113 mas (1.25 x 0.75 kpc) resolution shows continuum emission 
in BRI 1335-0417 with a total flux density of 208 ± 46 fiJy, consistent with 
the flux density measured with the VLA. The size of the source at FWHM is 
255 x 138 mas (1.7 x 0.9 kpc) and the derived intrinsic brightness temperature 
is ~ 3.5 x 10 4 K. No continuum emission is detected at the full VLBI resolution 
(32 x 7 mas, 211 x 46 pc), with a 4a point source upper limit of 34 fiJy beam -1 , 
or an upper limit to the intrinsic brightness temperature of 5.6 x 10 5 K. The 
highest angular resolution with at least a 4.5<r detection of the radio continuum 
emission is 53 x 27 mas (0.35 x 0.18 kpc). At this resolution, the image shows a 
continuum feature in BRI 1335-0417 with a size of 64 x 35 mas (0.42 x 0.23 kpc) 
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at FWHM, and intrinsic brightness temperature of ~ 2 x 10 5 K. The extent of the 
observed continuum sources at 1.4 GHz and the derived brightness temperatures 
show that the radio emission (and thus presumably the far-infrared emission) in 
BRI 1335-0417 is powered by a major starburst, with a massive star formation 
rate of order a few thousand M yr _1 . Moreover, the absence of any compact 
high-brightness temperature source suggests that there is no radio-loud AGN in 
this z = 4.4 QSO. 

Subject headings: galaxies: individual (BRI 1335-0417) — galaxies: active - 
galaxies: high-redshift — radio continuum: galaxies — techniques: interferomet- 
ric 



1. INTRODUCTION 



Optical surveys such as the Sl oan Digital Sky Survey (SDSS; York et al. 2000) and 
the Digitized Palomar Sky Survey ( Djorgovski et al. 19991) hay e re vealed large samples of 



Fan et al. 



(12002( 1 and iFan et all (12009 1 have 



quasi-stellar objects out to z ~ 6. Studies by 
shown that at such a high redshift we are approaching the epoch of reionization, the edge of 
the "dark ages" , when the first stars and massive black holes were formed. 

Observations of high redshift QSOs at mm and sub-mm wavelengths have shown that 
a significant fraction (~ 30%) of the sources are strong emitters of far-infrared (FIR) 
radiation, which is thermal emission from warm dust. These sources have luminosities 
L ftr > 10 12 L<ri, and mole cular gas masses greater than 10 10 M ( iSolomon fc Vanden Bout 



20051 : iRiechers et al.ll2006l ). Moreover, the large reservoirs of warm gas and dust in these ob 
jects have led to the hypothesis that th ese are starburst gal a xies with massive s t ar formation 
rates on the order of 10 00 M & yr" 1 faertoldi et aD 120031 : iBeelen et all bood : ICarilli et al" 
2004! : I Walter et al.ll2003h . 



An important question regarding these high- z QSOs is whether the dust heating mech- 
anism is dominated by a central AGN, or starburst activity in the QSO host galaxies. The 
high resolution of Very Long Baseline Interferometry (VLBI) observations permits a detailed 
look at the physical structures in the most distant cosmic sources. Also, sensitive VLBI con- 
tinuum observations can be used to determine the nature of the energy source (s) in these 
galaxies at radio frequencies. 

To date, several high redshift QSOs have been imaged at milliarcsecond re solution 
flFrev et alJ[l997l . l2003l : [Momiian et aDbood : lieelen et alJliooilMomiian et alihoosh . These 
are the highest resolution studies of such distant QSOs by far. In this paper, we present 
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sensitive VLBI observations of the z = 4.4 QSO BRI 1335-0417. 



The source BRI 1335-0417, which has an opti cal redshift of 4.396 ±0.026, was identified 
in th e Automatic Plate Measuring (APM) survey (jlrwin et al.lll99ll ; IStorrie-Lombardi et al. 
19961 ). This QSO is the second brightest source at millimeter wavele ngths with redshift s 
z > 4, with the brightest being the z = 4.7 QSO BRI 1202-0725 JOmont et al.lll996r i. 
The implied FIR luminosity of BRI 1335-0417 is L FIR = 3.1 x 10 13 L . Thus far, there 
i s no evidence for multiple imaging in BRI 1335-0417 due to strong gravitational lensing 
JCarilli et alill999h . 



The optical spectrum of the QSO BRI 1335-0417 shows a s trong absorption by Lya and 
metal lines at the source redshift (IStorrie-Lombardi et al.lll996f) . Its CO (5 — 4) line emission 
was detected at z = 4.4074±Q . 0015 (IGuilloteau et al.l 



1997J)- The CO (2-1) emission studies 



on this source by lCarilli et al.l (Il999l ) suggested optically thick emission from warm (> 30 K) 
molecular gas. The derived total molecular gas mass is Mm = (1.5 ±0.3) x 10 11 M , making 
BRI 1335-0417 one of the most massive molecular gas reservoirs known at high redshifts. 

The source BRI 1335-0417 is also detected in the radio continuum at 1.4 and 4.9 GHz 
with the Very Large Arr ay (VLA) in C-configuration, wit h flux densities of 220 ± 43 and 
76 ±11 yuJy, respectively (jCarilli et al.lll999l ; lYun et al.ll2000l ). The source appears unresolved 
at the resolution of these VLA observations, which is on the order of several arcseconds. 

Throughout this paper, we assume a flat cosmological model with fl m = 0.3, Qa = 0.7, 
and Hq = 71 km s -1 Mpc -1 . In this model, at the distance of BRI 1335-0417, 1 milliarcsec- 
ond (mas) corresponds to 6.6 pc. 



2. OBSERVATIONS AND DATA REDUCTION 

The VLBI observations of BRI 1335-0417 were carried out at 1.4 GHz on 2005 December 
31 and 2006 January 7, using the High Sensitivity Array (HSA), which included the Very 
Long Baseline Array (VLBA), the phased Very Large Array (VLA), and the Green Bank 
Telescope (GBT) of the NRAO0. Four adjacent 8 MHz baseband channel pairs were used in 
the observations, both with right and left-hand circular polarizations, and sampled at two 
bits. The data were correlated at the VLBA correlator in Socorro, NM, with 2 s correlator 
integration time. The total observing time was 14 hr. In these observations, the shortest 
baseline is between the phased VLA and the VLBA antenna in Pie Town, NM, which is 



1 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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52 km. This short-spacing limit filters out all spatial structure larger than about 0.42. 
Table 1 summarizes the parameters of these observations. 

The observations employed nodding-style phase referencing, using the calibrator J1335- 
0511 (S\_4 ghz — 0.4 Jy), with a cycle time of 4 min, 3 min on the target source and 1 min 
on the calibrator. A number of test cycles were also included to monitor the coherence of 
the phase referencing. These tests involved switching between two calibrators, the phase 
calibrator J1335-0511 and the phase-check calibrator J1332-0509 (Si. 4 gh z = 0.3 Jy), using 
a similar cycle time to that used for the target source. 

The accuracy of the phase calibrator position is important in phase-referencing obser- 
vations JWalkerill999l ). as this determines the accuracy of the absolute position of the target 



source and any associated components. Phase referen cing, as used her e, is known to preserve 



absolute astrometric positions to better than ±0. 01 ( Fomalontl 1 1 9991 ) . 



Data reduction and analysis were performed using the Astronomical Image Processing 
System (AIPS) and Astronomical Information Processing System (AIPS++) of the NRAO. 
After applying a priori flagging, amplitude calibration was performed using measurements 
of the antenna gain and system temperature for each station. Ionospheric corrections were 
applied using the AIPS task "TECOR". The phase calibrator J1335-0511 was self-calibrated 
in both phase and amplitude and imaged in an iterative cycle. 

Images of the phase-check calibrator, J1332-0509, were deconvolved using two different 
approaches: (a) by applying the phase and the amplitude self-calibration solutions of the 
phase reference source J1335-0511 (Figure la), and (b) by self calibrating J1332-0509 itself, 
in both phase and amplitude (Figure lb). The peak surface brightness ratio of the final 
images from the two approaches gives a measure of the effect of residual phase errors after 
phase referencing (i.e. 'the coherence' due to phase referencing). At all times, the coherence 
was found to be better than 98%. 

The self-calibration solutions of the phase calibrator, J1335-0511, were applied on the 
target source, BRI 1335-0417, which was then deconvolved and imaged at various spatial 
resolutions by tapering the visibility data. 



3. RESULTS & ANALYSIS 

Imaging the target source at the full resolution of the VLBI array, which is 32 x 7 mas 
(211 x 46 pc, PA=— 7°), achieved an rms noise level of 8.5 /iJy beam -1 , but did not reveal 
any continuum component in the field of BRI 1335-0417. This indicates the absence of any 
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compact radio continuum emission with flux densities of > 4a ~ 34 /ijy beam -1 , which in 
turn implies an upper limit to the intrinsic brightness temperature (corresponding to a rest 
frequency of ~8 GHz) of 5.6 x 10 5 K for any compact radio source in BRI 1335-0417. Our 
coherence tests during these observations using two VLBI calibrators show that the lack of 
a strong point source in BRI 1335-0417 at the full resolution of the array cannot be due 
to the phase referencing procedure. The 4cr flux limit reported above is almos t an order of 



magn itude lower than the flux measured by the VLA (220 ±43 /iJy) at 1.4 GHz (ICarilli et al. 



19991 ). This immediately implies that more than ~80% of the radio continuum emission in 



BRI 1335-0417 is extended and not confined to the central AGN. 

In the following, we will assess if our HSA observations can recover the flux seen by 
the VLA. To do so, we applied a two-dimensional Gaussian taper falling to 30% at 1 MA to 
the visibility data in both the u- and ^directions. This gives a beam size of 189 x 113 mas 
(1.25 x 0.75 kpc, P. A.= —48°), and the resulting image is shown in Figure 2. The rms noise 
level in this naturally weighted image is 28 /xJy beam -1 . The peak flux density of the detected 
continuum source is 126 ± 28 /xJy beam -1 , and the total flux density is 208 ± 46 //Jy, whic h 



agrees well with the flux density measured with the VLA at 1.4 GHz ( ICarilli et al.lll999l ). 
The size of the source at full width half maximum (FWHM) is 255 x 138 mas (1.7 x 0.9 kpc), 
and the derived intrinsic brightness temperature is (3.6 ± 0.9) x 10 4 K. 

Figure 3 shows the continuum emission in BRI 1335-0417 at the highest angular reso- 
lutions for which there is at least 4.5cr detection (a = 19 /iJy beam -1 ). This corresponds to 
a resolution of 53 x 27 mas (0.35 x 0.18 kpc) in position angle —6°. This image was obtained 
by applying a two-dimensional Gaussian taper falling to 30% at 5 MA in both the u- and 
^directions of the visibility data. The peak flux density of the continuum source detected 
at this resolution is 87 ± 19 /^Jy beam -1 , and the total flux density is 131 ± 30 /iJy. The 
size of the source is 64 x 35 mas (0.42 x 0.23 kpc) at FWHM, and its intrinsic brightness 
temperature is (2.2 ± 0.5) x 10 5 K. The results obtained at this resolution imply that about 
two thirds of the total radio emission emerges from the central ~0.3 kpc of BRI 1335-0417, 
and, likewise, one third arises from more extended scales (~1.3 kpc). 



4. DISCUSSION 



We have detected 1.4 GHz emission from BRI 1335-00417 (corresponding to a rest 
frame frequency of ~8 GHz) using the HSA. At a moderate resolution (189 x 113 mas; 
Figure 2) the intrinsic brightness temperature value of the detected continuum structure is 
~ 3.5 x 10 4 K. The me asured flux density at this resolution is consistent with the VLA flux 
density of this source (ICarilli et al.lll999l ). This implies that the radio continuum emission 
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at 1.4 GHz is confined to the extent of the structure seen in the VLBI image (Figure 2), 
which is 1.7 x 0.9 kpc at FWHM. Furthermore, the highest resolution image with continuum 
detection (Figure 3) shows that almost two thirds of the total radio continuum emission at 
1.4 GHz is arising from the central ~0.3 kpc. 

At the full resolution of our array (32 x 7 mas), the radio emission from BRI 1335- 
0417 is resolved out and does not show any single dominant source of very high brightness 



tempe rature (< 5.6 x 10 5 K). This is in contrast to the results obtained by iMomjian et al. 



(120041 ) on a sample of three high-z radio-loud quasars which were imaged with the VLBA, 
namely J1053-0016 (z = 4.29), J1235-0003 (z = 4.69), and J0913+5919 (z = 5.11). In each 
of these z > 4 quasars, a radio-loud AGN dominates the emission at 1.4 GHz on a few mas 
size scale, with intrinsic brightness temperatures in excess of 10 9 K. 



Condon et al.l (ll99ll ) have derived an empirical upper limit to the brightness tempera- 
ture for nuclear starbursts of ~ 10 5 K at 8 GHz (our rest frequency), while typical radio- loud 
AGN have brightness temperatures exceeding this value by at least two orders of magnitude. 
These authors also present a possible physical model for this limit involving a mixed non- 
thermal and thermal radio emitting (and absorbing) plasma, constrained by the radio-to-FIR 
correlation for star-forming galaxies. 

For BRI 1335-0417, the derived intrinsic brightness temperatures from our VLBI ob- 
servations are typical of starburst galaxies. Also, the linear extent of the radio continuum 
emission region in BRI 1335-0417 is typical of local starburst powered Ultra-Luminous IR 
galaxi es (ULIRGs; Sanders fc Mirabel 1996), such as Arp 220 Mrk 273, and IRAS 17208- 
0014 JSmith et alJll998l : ICarilli fc Tavloribood : iMomiian et al i 120031 ) . 



Another argument in favor of a massive starburst origin for the radio continuum emission 
from BRI 1335-0417 is that this source follows the radio-FIR correlation for star forming 
galaxies, although the mas sive star formation rate is extreme, of order a few thousand 
M yr" 1 JCarilli et alill999h . 



The physical characteristics of the QSO BRI 133 5-0417 are similar to two other high-z 
luminous IR sources; the brightest mm source in the lOmont et al.l (120031 ) sample of QSOs 
with redshifts be tween 1.8 < z < 2. 8, namely J1409+5628 at z = 2.58, and the brightest 
mm source in the lOmont et al.l (119961 ) sample of QSOs at redshifts z > 4, namely BRI 1202- 
0725 at z = 4.7. High angular resolution obser vations of J1409+5628 and BRI 1202-0725 at 
1.4 GHz (IBeelen et al.ll2004l ; ICarilli et al.ll2002l ) showed that the radio continuum emission in 
these optically very bright QSOs is dominated by extreme nuclear starbursts, with massive 
star formation rates on the order of 10 3 yr" 1 . M oreover, sensitive, tapered VLBI images 



of BRI 1202-0725 at 1.4 GHz (IMomjian et al.ll2005l ) revealed the starburst regions in each 
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of its two components that are separated by 4", with sizes and brightness temperatures 
comparable to that seen in BRI 1335-0417. 

The measured 1.4 GHz radio flux density in BRI 1335-0417 is consistent with extreme 
nuclear starburst activity, with the radio continuum being the sum of supernovae (SNe), 
supernova remnants, and re sidual relativist i c elec trons in the interstellar medium, as shown 
in the model pres ented bv ICondon et al.l (119911 ). However, detecting individual SNe at 
z = 4.4 is unlikely. ISmith et al.l (119981 ) reported the detection of several luminous radio SNe 
in the prototype ULIRG Arp 220 with flux densities between 0.2 and 1.17 mJy and angular 
sizes on the order of 0.25 mas. At the distance of BRI 1335-0417, the flux densities of such 
luminous radio SNe would be between (0.8 — 4.5) x 10 -3 fiJy. These values are three to four 
orders of magnitude lower that the rms noise levels achieved in our VLBI observations. 
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Fig. 1. — Continuum images of the phase-check calibrator J1332-0509 at 1.4 GHz: a) ob- 
tained by applying the phase and the amplitude self-calibration solutions of the phase refer- 
ence source J1335-0511, b) obtained by self calibrating J1332-0509 itself, in both phase and 
amplitude. The restoring beam size in both images is 27.9 x 6.5 mas in position angle —5°. 
The contour levels are at —3, 3, 6, 12, 24, 48 times the rms noise level in the phase-referenced 
image (top), which is 3.2 mJy beam -1 . The gray-scale range is indicated by the step wedge 
at the right side of each image. 
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Relative R. A. (mas) 

Fig. 2. — Naturally weighted 1.4 GHz continuum image of BRI 1335-0417 at 189 x 113 mas 
resolution (P. A. =—48°). The peak flux density is 126 fiJy beam" 1 , and the contour levels 
are at —2, 2, 3, 4 times the rms noise level, which is 28 //Jy beam" 1 . The gray-scale range 
is indicated by the step wedge at the right side of the image. The reference point (0, 0) 
is a(J2000.0) = 13 h 38 m 03?4034, 5(J2000.0) = -04°32'35'.'271. A two dimensional Gaussian 
taper falling to 30% at 1 MA was applied. 
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Fig. 3. — Naturally weighted 1.4 GHz continuum image of BRI 1335-0417 at 53 x 27 mas 
resolution (P. A. =—6°). The peak flux density is 87 ytxJy beam -1 , and the contour levels 
are at —2, 2, 3, 4 times the rms noise level, which is 19 yuJy beam -1 . The gray-scale range 
is indicated by the step wedge at the right side of the image. The reference point (0, 0) 
is a(J2000.0) = 13 h 38 m 03?4034, <5(J2000.0) = -04°32'35"271. A two dimensional Gaussian 
taper falling to 30% at 5 MA was applied. 
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Table 1. PARAMETERS OF THE VLBI OBSERVATIONS OF BRI 1335-0417 



Parameters Values 

Observing Dates 2005 Dec. 31 & 2006 Jan. 7 

Total observing time (hr) 14 

Phase calibrator J1335— 0511 

Phase-referencing cycle time (min) 4 

Frequency (GHz) 1.4 

Total bandwidth (MHz) 32 



